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We experimentally study the scattering of guided matter waves on an amplitude-modulated opti- 
cal lattice. We observe different types of frequency-dependent dips in the asymptotic output density 
distribution. Their positions are compared quantitatively with numerical simulations. A semiclassi- 
cal model that combines local Floquet-Bloch bands analysis and Landau-Zener transitions provides 
a simple picture of the observed phenomena in terms of elementary Floquet photon absorption- 
emission processes and envelope-induced reflections. Finally, we propose and demonstrate the use 
of this technique with a bichromatic modulation to design a tunable sub-recoil velocity filter. Such 
a filter can be transposed to all species since it does not rely on a specific internal level configuration 
of the atoms. 

PACS numbers: 03.75.Kk,03.75.Lm 



I. INTRODUCTION 

Cold atoms interacting with time-modulated optical 
lattices display a wide variety of quantum and classical 
dynamics. These include the observations of dynami- 
cal localization [TJ [2] , chaos- assisted tunneling [3, 4 , the 
Anderson metal-insulator transition in momentum space 
[5], dynamically controlled tunneling [6HH] , matter wave 
engineering [9HT2] or the probing of gapped modes in de- 
generate quantum gases JTSJ [14] . 

In this article we propose and demonstrate the use of 
time-dependent optical lattice for atom optics through 
the investigation of the scattering of a cold atom packet 
on an amplitude-modulated optical lattice [15] . Our tech- 
nique provides a new tunable tool for velocity selection in 
confined geometries and a system for studying quantum 
transport with time-dependent potentials [16] as initially 
studied in mesoscopic physics [T71 [18] . 

The Bragg reflection of a propagating guided matter 
wave on a static optical lattice has been recently demon- 
strated [19] . The Bragg condition on the class of veloci- 
ties that are reflected reads ^Bragg = where n is an 
integer, v\, = h/(md) and d is the lattice spacing. This 
condition is valid in the perturbative limit i.e. for a small- 
depth optical lattice Uo <C El where El = mv L /2 is the 
lattice energy scale. However ^Bragg is directly related 
to the periodicity of the lattice and therefore cannot be 
easily tuned over a large range. 

In contrast, the interaction of a propagating matter 
wave with an amplitude-modulated optical lattice of- 
fers more flexibility and realizes, in particular, a tunable 
Bragg reflector. Our study explores the non perturbative 
regime for which the lattice depth is on the order of the 
lattice energy scale El and the modulation depth is rela- 
tively large. In this range of parameters and for our typ- 
ical interaction time, the stroboscopic phase space does 
not exhibit any chaos [20 and the scattering can there- 
fore be used to engineer the velocity distribution of the 
incoming packet of atoms. 

The paper is organized as follows. In Sees. [TT] and |III[ 



we present the experimental setup and review the main 
observations. In Sec. |TV} we derive a simple model based 
on a modulated vanishing depth optical lattice. This ap- 
proach is successfully compared with the experimental 
results. In Sees. W\ and [VT] we develop a semiclassical 



model that combines the local Floquet-Bloch framework 
and Landau-Zener transitions to get a deeper insight on 
the elementary processes at work in the matter wave scat- 
tering. In Sec. |VII[ we present an application with the 
realization of a tunable notch velocity filter through the 
scattering on an optical lattice whose amplitude is time- 
modulated by a bichromatic field. Conclusions are drawn 

in Sec. gnu 



II. EXPERIMENTAL SETUP 

Our experimental setup has already been described 
in [19]. In short, a thermal cloud of typically 8.10 4 
rubidium-87 atoms at T — 500 nK is obtained after 3.5 s 
of forced evaporation in a crossed dipole trap formed by 
two red-detuned (1070 nm) laser beams: a horizontal 
guide and a dimple beam. During the evaporation, we 
use the spin-distillation technique to prepare atoms in 
\F = l,m F = 0) [2U [22]. We deliberately use a thermal 
cloud rather than a BEC in order to probe the mod- 
ulated lattice for a wide range of velocities in a single 
shot (see below). By switching off the dimple beam, 
we release a packet of longitudinal velocity dispersion 
Av = 6 mm/s in the horizontal guide. Atoms are sub- 
sequently accelerated by a t acc = 15 ms inhomogeneous 
magnetic field pulse to a mean velocity v = 10 mm/s. 
The atomic packet then propagates towards the lattice 
whose center is located 500 fim downstream from the 
trap position (see Fig. [TJa)). The lattice is obtained by 
crossing two horizontal off-resonance laser beams (wave- 
length A = 840 nm, waist w = 100 //m) at an angle 
= 81° [19]. We modulate the lattice intensity using an 
acousto-optic modulator prior to the beam separation. 
The time-dependent potential experienced by the atoms 
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FIG. 1: (color online) (a) Sketch of a propagating atomic 
packet impinging onto an optical lattice whose amplitude is 
modulated, (b) Measured longitudinal density n(x,t acc + 
tprop) after an acceleration stage of t acc = 15 ms and a prop- 
agation time tprop = 78 ms for various lattice modulation 
frequencies v (lattice depth Uo = 2El, modulation depth 
rj = 33%, lattice position given by the dotted line), (c) Nu- 
merical simulations (see text) with a resolution that matches 
the experimental optical resolution (~ 10 /im). Frequency- 
dependent dips are observed in the transmitted distribution. 
Dotted, dot-dashed and dashed lines in (c) show the linear 
dependence of the dip position with v. 



reads: 

U(x^ = -Uo(t)e- 2x2 ^\m"{^), (1) 

where U (t) = U (1 + rj cos (2nut)), d = A/[2sin(0/2)] = 
650 ± 15 nm (v^ = 7.1 mm/s, El/H — vl = 5.4 kHz) 
and w = w/ cos(#/2) ~ 130 /im. The lattice depth Uo = 
2El is calibrated by Kapitza-Dirac diffraction [23 . The 
typical modulation depth is r] = 33%. The atomic packet 
propagates during t prop = 78 ms through the lattice and 
is imaged in-situ (without time-of- flight) by absorption 
imaging. 



III. EXPERIMENTAL RESULTS 

Figure [T](b) shows the atomic density after propagation 
as a function of the modulation frequency v. Each hori- 
zontal line is obtained by averaging 8 images integrated 
along the transverse direction. Two kinds of density dips 



are observed in the transmitted part: (i) dips whose po- 
sitions do not depend on the modulation frequency and 
that correspond to velocity classes fulfilling the Bragg 
reflection on the static lattice [19 , and (ii) dips whose 
positions depend on the frequency. As we shall discuss 
below, some dips of the latter category have their coun- 
terpart in the reflected packet and correspond to reflected 
class of velocity while others are due to slowing down or 
acceleration effects. 

Except for the zones very close to depletion lines in the 
transmitted part in Fig. [TJb), each position downward 
the lattice can be mapped onto a well denned class of 
incident velocity x c± Vi nc (^acc + ^prop) + K where K is a 
constant [24 . Let us characterize the different depletion 
lines. Using the correspondence between x and v- inci the 
main depletion line (white dashed line) in Fig. [TJc) has 
a slope (660 nm) -1 ~ 1/d. We also observe directly 
the corresponding reflected atoms in the region x < 0. 
The red upper dashed line of depleted atoms in Fig. [TJc) 
is parallel to the main line and has no counterpart in 
the reflected region. The white dot-dashed (dotted) line 
in Fig. [TJc) has a slope twice (three times) as large as 
the one of the dashed white and red depletion lines in 
Fig. 0c). 

Figure [TJc) is the result of a numerical simulation of 
the atomic packet dynamics using the one-dimensional 
Schrodinger equation solved by the split-Fourier method 
and a wave packet whose initial momentum and position 
dispersions match the measured experimental values [25J. 
We find a good agreement between simulations and ex- 
periment. However such an approach does not reveal the 
mechanisms responsible for the depletion lines. 



IV. THE MODULATED VANISHING DEPTH 
MODEL 

The slopes and the relative position of the depletion 
lines can be simply interpreted in terms of interband 
transitions in the limit of a small lattice depth. Consider 
an incident quasi-monochromatic wave packet of velocity 
vq. This velocity dictates the band in which the incoming 
atom enters. For instance, if 3^l/2 < vq < 2v\, the atom 
will be on band 4 (see Fig. |2|. 

In the limit of a vanishing lattice depth, the band 
structure can be constructed by the superposition of 
parabolic energy spectra centered around all reciprocal 
points E n (k) = h 2 (k — n/cz,) 2 /2m, where n is an inte- 
ger and kL = 2ir/d. For instance, bands 2 and 3 are 
constructed from the parabolas centered at ±hkL- 

In the limit of a small amplitude modulation, the mod- 
ulation drives interband transitions that keep the pseudo- 
momentum hk unchanged. In Fig.[3j we have represented 
schematically the experimentally observed depletion lines 
in velocity space (see Fig.JT]). For each depletion lines and 
in the different velocity domains, we have also indicated 
two numbers: the first corresponds to the band on which 
atoms lies according to their incident energy E = ttivq/2 
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FIG. 3: (color online) Schematics of the different depletion 
lines observed in Fig. 1 of the article along with the numbers 
of the bands involved in the transition dirven by the frequency 
of the amplitude modulation. With the experimental param- 
eters, — 7.1 mm/s. 



FIG. 2: The band diagram is constructed by the superpo- 
sition of parabolic energy spectra centered around all recip- 
rocal points E n (k) = h 2 (k — n/cz,) 2 /2ra. For each band we 
have indicated the index n of the parabola from which the 
corresponding branches have been constructed. 

while the second corresponds to the band on which atoms 
are resonantly transferred by the modulation. 

The interband transition frequencies that promote an 
atom from parabola n to n' are given by 

±i/ n _> n / = (E n -E n ,)/h. (2) 

The sign + (— ) corresponds to a transition to a lower 
(upper) band. By energy conservation, we have 

^mvl = E u (k). (3) 

Combining Eqs. Q and we get 

E n -E n > /2 n'-n 
± v n ^ n > = = -{n -n) vl -\ — v . (4) 

As a first example, let us consider the green dotted 
depletion line (see Fig. [3|. In the incident velocity do- 
main 3^l/2 < < 2^l, it corresponds to a resonant 
transition between bands 4 and 3. According to Fig. [2] 
we have a priori two possibilities either a transition from 
parabola n = 2 to parabola n' = — 1 or a transition from 
parabola n = — 2 to parabola n r = 1 occurring at ±ac. 



The slope of branch n — — 2 is positive which thus corre- 
sponds to a propagation from left to right (positive group 
velocity). The states of the incoming packet with posi- 
tive velocity in the range of energy that corresponds to 
band 4 will thus be projected onto the states that cor- 
responds to the branch n = — 2. Knowing n and n' we 
deduce from Eq. Q, the equation for the dotted line is 
v — —9vl + 3vo/d. 

Such an interpretation can be made for all depletion 
lines. The lower (black) dashed line of Fig. [3] corresponds 
successively to the interband transitions 3 — >• 1, 4 — >• 2, 
5 — >• 3 ... From Eq. Q, we deduce its equation v = 
— vl -\-vo/d. Similarly we find for the upper (red) dashed 
line v = vjj + vq/oI. The offset between the two dashed 
lines is thus equal to 2v^ ~ 10.8 kHz. We obtain for the 
(blue) dot-dashed line the equation v = —Avl + 2vo/d. 

To get a better understanding of the width of the de- 
pletion lines, their interpretation in terms of elementary 
processes, the timescale on which the transitions occur 
and the role played by the Gaussian envelope of the lat- 
tice potential, we introduce now a more elaborated anal- 
ysis based on Floquet-Bloch framework [26l [27] . 



V. THE FLOQUET-BLOCH FRAMEWORK 

This approach is not restricted to small modulation 
depth and is thus well-adapted to analyze the experimen- 
tal situation. For a potential periodic in both space and 
time, the Floquet-Bloch solutions of the time-dependent 
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FIG. 4: (color online) (a) Floquet-Bloch band diagram for 
a square- envelope lattice (depth Uo = 2El) and modulation 
frequency v — 11 kHz. Band color code: green —0.5 < (uf) < 
0.5, red (tif) > 0.5 and blue (uf) < —0.5 where uf is the 
Floquet excitation number, (b) Probability of reflection ob- 
tained from a numerical simulation of a ID wavepacket with 
an incident velocity dispersion Av — 0.2 mm/s impinging 
on time amplitude-modulated lattice with a finite square- 
envelope (length=80 d, Uo = 2El, rj = 30%) as a function 
of the incident energy and v. The horizontal white dashed 
line shows the case v — 11 kHz that corresponds to the dia- 
gram (a). The horizontal dashed (dotted) lines denotes open 
gap (degenerate) anticrossings. Only open gap anticrossings 
yield reflection. 



Schrodinger equation read: 

t/> n , k (x, t) = e^-^l^un^x, t), (5) 

where e n (k) are the quasi-energies. The functions 
u ny k(x,t) are biperiodic in space and time and therefore 
can be Fourier expanded: 



t1 k(x,t) = u nyk (x + d,t) 



EE< 



k(x,t + T) 

F ut) 



(6) 



In the following, we restrict ourselves to rip G { — 1,0,1} 
i.e. to situations in which only one Floquet photon can be 
absorbed or emitted [28]. At zero modulation depth, the 
Floquet-Bloch band diagram is nothing but the superpo- 
sition of the Bloch diagrams shifted by nphuj. At finite 
modulation depth, anticrossings appear for frequencies 
that correspond to interband transitions. 

Consider first the simple case of a square-envelope lat- 
tice of amplitude Uo modulated at a frequency v with 
an amplitude rj. In Fig. [4ja) we plot the Floquet-Bloch 
spectrum for v = 11 kHz. Two kinds of anticrossings can 
be identified: those yielding open gaps (horizontal dashed 
line in Fig.Q and those without gaps for which two states 
with the same quasi-energy are available (horizontal dot- 
ted line in Fig. [4|. To identify the role of the different 
types of anticrossings on the incident matter wave packet, 
we have performed a ID simulation that solves the corre- 
sponding time-dependent Schrodinger equation. Figure 
[4^b) gives the reflection coefficient as a function of the in- 
cident energy Eq and the modulation frequency v. Two 
types of reflection conditions can be clearly identified: (i) 



those due to Bragg reflection onto the static lattice (no 
dependence on v) and (ii) those that correspond to open 
gap anticrossings and whose positions depend on v. The 
interpretation is clear: when the incident energy falls in 
an open gap anticrossing, no propagating state is avail- 
able and the particle is reflected. The degenerate anti- 
crossings do not induce reflection in the square-envelope 
case. However, as we discuss below, they play an impor- 
tant role in the dynamics of the experimentally relevant 
case in which the lattice has a slowly varying envelope. 



VI. THE SEMICLASSICAL MODEL 

In this latter case, the situation turns out to be dif- 
ferent since the system can follow adiabatically a quasi- 
energy band during its time evolution. To describe this 
propagation, we propose a semiclassical model that en- 
ables one to identify the elementary processes responsi- 
bles for the velocity changes of the particle and the posi- 
tion at which such processes occur. It contains two main 
ingredients. The first one consists in describing the par- 
ticle motion on a given local Floquet-Bloch band through 
the combined evolution of the wavepacket position and 
of its mean pseudo momentum k. The corresponding set 
of coupled equations reads 



lcten 
h dk 



and 



k 



1 de n 
h dx ' 



(7) 



The first equation defines the group velocity of the wave 
packets while the second results from the adiabatic fol- 
lowing condition de n (k,x)/dt = [26j[29]. The second 
ingredient consists in taking into account the possibility 
for a particle to undergo a Landau-Zener transition when 
it passes an anticrossing. In our case, the approximation 
of a local two-level situation is valid and therefore the 
probability to change the band index is P = e _27r7 with 



7 : 



AE 2 

Ah 



dt 



(e n — £ n ±i) 



(8) 



where AE is the size of the gap [30 , 3T]. 

In practice, we evolve the particle according to Eqs. ^ 
and we evaluate the energy difference to lower or upper 
bands at each time step. When this quantity reaches a 
minimum (i.e. at an avoided crossing position), we com- 
pute the corresponding Landau-Zener probability and 
transfer or not the particle to the next band according 
to this probability. 

To validate this semiclassical trajectory method, we 
compare it with the full resolution of the corresponding 
ID Schrodinger equation. To perform this comparison, 
we have simulated the semiclassical trajectories of 1700 
incoming velocities about the mean velocity of the packet 
in the following range —4.5 mm/s < V[ nc < 24.5 mm/s. 
In this way, we sample 98 % of the initial distribution. 
Furthermore, we perform 35 shots for each incoming ve- 
locity to improve the statistics of our Monte Carlo simu- 
lation. 
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FIG. 5: Dashed line: Result of the resolution of the ID 
Schodinger equation. Solid line, result of the semiclassical 
simulation with random Landau-Zener transition. The ar- 
rows indicate the main processes for three depletion lines. 
The wave packet has a mean velocity v — 10 mm/s and a 
velocity dispersion Av — 6 mm/s. 



For each velocity class v- inci we get the density from 
the final positions of the different shots weighted by the 
initial wavepacket velocity distibution density at ^i nc . 
Figure [5] provides an example of such a comparison for 
v = 20 kHz. The key features are very well captured by 
the semiclassical simulation. In this plot, we can easily 
identify the depletion lines that correspond to the differ- 
ent interband transitions. 

We can use the semiclassical model to analyze the 
different mechanisms yielding to depletion bands as ob- 
served in Fig. [I] For this purpose, we follow determinis- 
tically the branch for which the Landau-Zener transition 
probability is above 1/2 at an avoided crossing. To il- 
lustrate the wide variety of possibilities, we shall choose 
three generic and different set of parameters v) yield- 
ing to dips in the output density distribution (see labels 
1, 2 and 3 in Fig. [ljc)). In Fig. |6ja), we plot the veloc- 
ity along with the mean Floquet excitation number for 
each case and for the main trajectory given by the Monte- 
Carlo simulation. In Fig.[6|b) we show the corresponding 
Floquet-Bloch diagrams in the region of interest. When 
a particle is moving toward the center, all quasi ener- 
gies decrease since the amplitude of the attractive lattice 
increases. As a result, the particle state moves up rel- 
atively to the band diagram. In the same way, if the 
particle is moving backward, it will go down the hills of 
the diagram. With these simple pictures in mind trajec- 
tories can be readily interpreted. In each case, the key 
phenomenon is the absorption or stimulated emission of 
a Floquet photon by adiabatic following denoted by a in 
Fig. [6] In case (1), the particle emits a Floquet photon, 
performs a reflection (denoted (3) when reaching the bot- 
tom of the band and absorbs a Floquet photon (a!) be- 
fore leaving the lattice. In case (2), the first emission only 



slows down the particle which is then Bragg reflected and 
subsequently accelerated by Floquet photon absorption. 
In case (3), the particle is not reflected. It is transiently 
accelerated in the lattice by a Floquet photon absorption- 
emission cycle. For a much longer propagation time, the 
dip would be refilled. 

Other features of the experimental and numerical di- 
agrams of Fig. [I] can be readily explained thanks to our 
semiclassical model. For instance, the density bump 
above the white dashed line corresponds to atoms that 
have been slowed down. The velocity spread of reflected 
particles at position 2 (see Fig. [I]) is 0.8 mm/s. This 
value can be recovered from the variation of the energy 
position of the gap along the lattice. 

For a given incident kinetic energy Eq, a large size of 
the envelope and/or a large modulation depth increases 
the efficiency of the process since it favors an adiabatic 
following of the anticrossings. A less intuitive feature 
concerns the lattice depth. Indeed, a small lattice depth 
(Uo < Eq) increases the selectivity of the class of incident 
velocities that are affected by the modulation. This origi- 
nates from the fact that the system is projected on a high 
energy band, and the position of the gap remains roughly 
constant throughout the lattice. Interestingly enough, 
this ensures the robustness of the method against the 
specific shape of a smooth envelope. 



VII. APPLICATION: VELOCITY FILTER 

The narrowest velocity filters used in the cold atom 
community rely on velocity selective Raman transitions 
on atoms in free space. This technique involves a com- 
bined change of internal and external states. The achiev- 
able velocity width are in the range of 200-300 fim/s 
IH EL 38J . Using the scattering on amplitude modulated 
optical lattice we demonstrate hereafter a new technique 
to realize a velocity filter with a width slightly larger than 
the state of the art with velocity selective Raman tran- 
sitions. Our technique uses only the external degrees of 
freedom and thus does not require any specific internal 
configuration. In addition it is well adapted for guided 
matter waves. 

We turn our device into a tunable momentum filter by 
combining different modulation frequencies. We use here 
the main reflection line (white dashed line in Fig. [I]) that 
acts as a notch filter in momentum space. For this pur- 
pose, we modulate the lattice with two different frequen- 
cies to create a transmitted band between two rejected 
ones: Uo(t) = Uo(l+r] cos(27rz/it)+?7 cos(27r^2^))- Strictly 
speaking, the detailed dynamics of a wave packet submit- 
ted to this two frequency and non perturbative modula- 
tion cannot be inferred directly from the single frequency 
dynamics [33 . However, the simple picture according to 
which nearly independent dips can be drilled into the ve- 
locity distribution with two frequencies is quite robust. 
We observe that the reflection spectrum is roughly the 
product of the two independent spectra (see Fig. [7| [34] . 
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FIG. 6: (color online) (a) Velocity and mean Floquet excitation number (uf) as a function of the propagating time for 
parameters corresponding to three different depletion zones shown in Fig. [I] (b) Local Floquet-Bloch diagram. Dark arrow 
denotes the trajectory followed by the fictitious particle of the semiclassical model (see text). Case (1): Reflection on an open 
gap v = 10.3 mm/s, v — 11 kHz. Case (2): Reflection on an anticrossing without gap v = 15.8 mm/s, v — 20 kHz. Case (3) 
transient acceleration v = 15.8 mm/s, v — 30.5 kHz. a (f3) denotes absorption or emission of one Floquet photon (reflection). 




FIG. 7: (color online) (a) Density distribution for the scatter- 
ing on a two-frequency modulated optical lattice (same con- 
ditions as in Fig. [TJ for a fixed frequency v\ — 16 kHz and 
a scanned frequency v^. A narrow slice of transmitted atoms 
is produced (see arrow), (b) Number of atoms and velocity 
dispersion associated with these narrow slices of transmitted 
atoms as a function of V2- Inset: mean velocity of the slice of 
atoms as a function of z/ 2 . 



velocity dispersion of the velocity filter that we have de- 
signed is on the order of 450 /im/s for our parameters 
(i.e. 1.1 nK in temperature units). 

VIII. CONCLUSION 

The matter wave engineering presented here does not 
have any fundamental limit. A further improvement of 
velocity selection could be achieved using a smaller depth 
lattice combined with a larger waist size. This technique 
can be transposed easily to other species since it does 
not rely on a specific internal level configuration. By con- 
struction, it is well adapted to ID geometry and therefore 
enhances the toolbox of guided atom optics. Finally, the 
control of the guided atomic flux for a given and tunable 
narrow class of velocities as studied here is reminiscent 
of the quantum modulated transistor principle where the 
gate voltage is replaced by the modulation [39 . 



The mean velocity of the slice of atoms is therefore gov- 
erned by d(y\ + 1*2) /2 while its width is controlled by the 
frequency difference \v<i — v\\. In our set of experiments, 
v\ is fixed at 16 kHz and V2 is varied from 16 to 20 kHz. 
Between the two reflection lines, atoms in a narrow class 
of velocity are transmitted (arrow in Fig. [7]). The slice 
contains about 1000 atoms and has a mean velocity on 
the order of 15 mm/s (inset of Fig. [7]). The minimum 



Acknowledgments 

We thank useful discussion with R. Dubertrand, G. 
Lemarie and B. Georgeot. We acknowledge financial 
support from the Agence Nationale pour la Recherche, 
the Region Midi-Pyrenees, the university Paul Sabatier 
(OMASYC project), the NEXT project ENCOQUAM 
and the Institut Univer sit aire de France. 



[1] B. G. Klappauf, W. H. Oskay, D. A. Steck, and M. G. [2] J. Ringot, P. Szriftgiser, J. C. Garreau, and D. Delande, 
Raizen, Phys. Rev. Lett. 81, 1203 (1998). 



7 



Phys. Rev. Lett. 85, 2741 (2000). 
[3] W. K. Hensinger et a/., Nature 412, 52 (2001). 
[4] D. A. Steck, W. H. Oskay and M. G. Raizen, Science 

293, 274 (2001). 
[5] J. Chabe, G. Lemarie, B. Gremaud, D. Delande, P. Szrift- 

giser, and J. C. Garreau, Phys. Rev. Lett. 101, 255702 

(2008). 

[6] H. Lignier, C. Sias, D. Ciampini, Y. Singh, A. Zen- 
esini, O. Morsch, and E. Arimondo, Phys. Rev. Lett. 99, 
220403 (2007). 

[7] E. Kierig, U. Schnorrberger, A. Schietinger, J. Tomkovic, 
and M. K. Oberthaler, Phys. Rev. Lett. 100, 190405 
(2008). 

[8] Y.-A. Chen, S. Nascimbene, M. Aidelsburger, M. Atala, 
S. Trotzky, and I. Bloch, Phys. Rev. Lett. 107, 210405 
(2011). 

[9] A. Alberti, G. Ferrari, V. V. Ivanov, M. L. Chiofalo and 
G. M. Tino, New J. Phys. 12, 065037 (2010) . 
[10] S. Arlinghaus and M. Holthaus, Phys. Rev. A 84, 063617 
(2011). 

[11] C. E. Creffield and F. Sols, Phys. Rev. A 84, 023630 
(2011). 

[12] Jacob F. Sherson, Sung Jong Park, Poul L. Pedersen, 

Nils Winter, Miroslav Gajdacz, Sune Mai and Jan Arlt, 

New J. Phys. 14, 083013 (2012). 
[13] R. Jordens, N. Strohmaier, K. Giinter, H. Moritz and T. 

Esslinger, Nature 455, 204 (2008). 
[14] S. D. Huber, B. Theiler, E. Altman, and G. Blatter, Phys. 

Rev. Lett. 100, 050404 (2008). 
[15] M. Henseler, T. Dittrich, and K. Richter, Phys. Rev. E 

64, 046218 (2001). 
[16] M. Biittiker and R. Landauer, Phys. Rev. Lett. 49, 1739 

(1982). 

[17] R. H. Blick, D. W. van der Weide, R. J. Haug, and K. 

Eberl, Phys. Rev. Lett. 81, 689 (1998). 
[18] G. Platero and R. Aguado, Phys. Rep. 395, 1 (2004). 
[19] C. M. Fabre, P. Cheiney, G. L. Gattobigio, F. Vermersch, 

S. Faure, R. Mathevet, T. Lahaye, and D. Guery-Odelin, 

Phys. Rev. Lett. 107, 230401 (2011). 
[20] Private communication from R. Dubertrand, G. Lemarie 

and B. Georgeot. 
[21] A. Couvert, M. Jeppesen, T. Kawalec, G. Reinaudi, 

R. Mathevet and D. Guery-Odelin, Europhys. Lett. 83, 



50001 (2008). 

[22] G. L. Gattobigio, A. Couvert, M. Jeppesen, R. Math- 
evet, and D. Guery-Odelin, Phys. Rev. A 80, 041605(R) 
(2009). 

[23] Yu. B. Ovchinnikov, J. H. Miiller, M. R. Doery, E. J. 

D. Vredenbregt, K. Helmerson, S. L. Rolston, and W. D. 

Phillips, Phys. Rev. Lett. 83, 284 (1999). 
[24] If a denotes the acceleration experienced by the atoms 

during the magnetic pulse, one finds K = at &cc (t &cc /2 — 

tpvog)- 

[25] In such a scattering experiment, atom- atom interactions 
are negligible in the course of the propagation. 

[26] N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(W. B. Saunders Co., Philadelphia, 1976). 

[27] Jon H. Shirley, Phys. Rev. 138, B979 (1965). 

[28] This assumption suffices to capture all observed phenom- 
ena on the experiment. The extension to larger Floquet 
excitations is straightforward. 

[29] K. Drese and M. Holthaus, Eur. Phys. J. D 5, 119 (1999). 

[30] L. D. Landau and E. M. Lifshitz, Quantum mechanics: 
Non-relativistic theory, 3rd ed. (Pergamon, New York, 
1977), pp. 342-351. 

[31] C. Zener, Proc. R. Soc. London Ser. A137, 696 (1932). 

[32] The asymptotic values that enter Eq. (|8| are obtained 
by a quadratic approximation of the energy difference 
\s n — £ n ±i\- This ensures the robustness of the calculation 
of the Landau Zener probability. 

[33] S.-I. Chu and D. A. Telnov, Physics Reports 390, 1 
(2004). 

[34] The simple picture of two independent frequencies breaks 

down when the frequency difference is below 0.5 kHz. 
[35] P. L. Gould, G. A. Ruff and D. E. Pritchard, Phys. Rev. 

Lett., 56, 827 (1986). 
[36] D. W. Keith, M. L. Schattenburg, Henry I. Smith, and 

D. E. Pritchard, Phys. Rev. Lett. 61, 1580 (1988). 
[37] M. Kasevich and S. Chu, Phys. Rev. Lett. 69, 1741 

(1992). 

[38] R. Battesti, P. Clade, S. Guellati-Khelifa, C. Schwob, B. 

Gremaud, F. Nez, L. Julien and F. Biraben, Phys. Rev. 

Lett. 92, 253001 (2004). 
[39] F. Sols, M. Macucci, U. Ravaioli and K. Hess, J. Appl. 

Phys. 66, 3892 (1989). 



